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Executive  Summary 


This  research  effort  represents  our  investigation  on  the  following  problems: 

1.  Optimal  processor  structure  for  detecting  a  rank  one  signal  in  additive  non-Gaussian 
radar  modeled  by  a  spherically  invariant  random  process  (SIRP)  and  its  performance 
analysis. 

2.  Analytical  expressions  for  the  detection  and  false  alarm  probability  for  two  impor¬ 
tant  detection  tests. 

3.  Statistical  analysis  of  the  non-homogeneity  detector  for  Gaussian  interference  sce¬ 
narios. 

Our  work  on  item  1  presents  integral  expressions  for  detection  and  false  alarm  prob¬ 
abilities  of  the  optimal  processor  for  detecting  a  rank  one  signal  in  a  background  of 
non-Gaussian  clutter  which  is  modeled  by  a  SIRP.  It  is  well  known  that  when  operating 
in  Gaussian  clutter  backgrounds,  the  best  achievable  performance  is  that  of  a  matched 
filter  (MF)  or  a  normalized  correlator.  In  this  work  we  present  the  performance  analysis 
of  the  corresponding  optimal  processor  for  SIRPs.  Integral  expressions  for  detection  and 
false  alarm  probabilities  are  obtained  [1]. 

In  many  instances,  it  is  impossible  to  implement  the  optimal  processor  outlined  in 
item  1.  Consequently,  we  seek  methods  which  circumvent  the  need  to  know  the  clutter 
probability  density  function  (PDF).  This  aspect  is  addressed  in  item  2,  where  we  present 
analytical  results  for  the  probability  of  detection  and  probability  of  false  alarm  for  the 
normalized  matched  filter  (NMF)  and  normalized  adaptive  matched  filter  (NAMF)  tests. 
Performance  is  compared  to  the  adaptive  matched  filter  and  the  Kelly  GLRT  receiver. 
The  analytic  expressions  derived  in  this  work  were  reported  in  [2]. 

Finally,  an  important  issue  in  STAP  is  that  of  homogeneity  of  training  data.  Non¬ 
homogeneity  of  the  training  data  has  a  deleterious  effect  on  STAP  performance  in  that 
undernulled  clutter  significantly  degrades  detection  and  false  alarm  characteristics.  Previ¬ 
ous  work  in  this  area  has  proposed  the  use  of  a  non-homogeneity  detector  (NHD)  based  on 
a  generalized  inner  product  (GIP).  Our  work  on  item  3  considers  the  statistical  analysis 
of  the  NHD  for  Gaussian  interference  statistics.  We  show  that  a  more  stringent  test  can 
be  constructed  by  accounting  for  the  statistics  of  the  (GIP )  under  the  condition  of  finite 
training  data  support.  In  particular,  exact  theoretical  expressions  for  the  GIP  PDF  and 
GIP  mean  are  derived.  Additionally,  we  show  that  for  Gaussian  interference  statistics, 
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based  test  is  presented.  Results  from  this  investigation  were  reported  m  [3]. 
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Chapter  1 


Optimal  Space-Time  Adaptive 
Processing  Method  in  Non-Gaussian 
Radar  Clutter  Backgrounds 


1.1  Motivation 

This  analysis  is  motivated  by  the  problem  of  space-time  adaptive  (STAP)  processing  m 
non-Gaussian  clutter  backgrounds.  The  problem  of  adaptive  target  detection  in  Gaussian 
clutter  has  received  considerable  attention  [4-8].  The  corresponding  prob  em  or  non 
Gaussian  clutter  backgrounds  has  been  the  focus  of  recent  work  reported  in  19-14].  ine 
receiver  developed  in  [7, 15]  is  an  adaptive  matched  filter  (AMF),  whose  performance 
asymptotically  reaches  that  of  the  matched  filter  (MF)  which  is  the  optimal  processor  for 
Gaussian  clutter  statistics.  Thus,  the  MF  represents  an  upper  bound  on  the  performance 
of  the  AMF  for  Gaussian  clutter  statistics. 

In  [10, 11],  we  showed  that  the  receiver  for  non-Gaussian  SIRPs  is  equivalent  to  an 
adaptive  matched  filter  compared  to  a  data  dependent  threshold.  The  distribution  of  the 
test  statistic  was  extremely  difficult  to  derive  in  closed  form.  Consequently,  performance 
analysis  in  [11]  was  carried  out  by  Monte-Carlo  simulation  in  the  multichannel  signal 
processing  system  developed  at  AFRL/SNRT.  The  goal  of  this  effort  is  to  obtain  an  upper 
bound  on  the  performance  of  the  adaptive  processor  provided  by  the  optimal  receiver  for 
SIRPs.  This  enables  the  specification  of  meaningful  adaptive  performance  metrics  such  as 
the  sample  support  size  needed  to  get  to  within  3dB  of  the  optimal  processor  performance. 

Accordingly,  we  briefly  outline  the  derivation  of  the  optimal  processor  for  SIRPs  and 
present  relevant  details  of  the  performance  evaluation.  Integral  expressions  for  the  proba¬ 
bility  of  detection  and  false  alarm  are  derived.  As  a  special  case,  these  expressions  reduce 
to  the  well  known  Gaussian  MF  results. 
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1.2  Problem  Statement 


We  consider  the  following  statistical  hypothesis  testing  problem: 

Ho  :  x  =  y  (1.1) 

Hi :  x  =  as  +  y 


where  x  =  Received  JN  x  1  complex  observation  vector 
s  =  Known  steering  vector 
a  =  Unknown  complex  signal  amplitude 

y  =  Complex  spherically  invariant  random  vector  (SIRV)  with  known  Hermitian  positive 
definite  covariance  matrix,  £  and  known  characteristic  probability  density  function  (PDF) 
fv(v).  It  is  important  to  note  that  the  optimal  processor  for  SIRVs  assumes  knowledge 
of  the  covariance  matrix  and  the  characteristic  PDF. 

The  maximum  likelihood  estimate  of  a  is  given  by  [11] 


.  sH£_1x 


(1.2) 


This  estimate  is  then  used  in  a  likelihood  ratio  test  to  discriminate  between  the  two 
hypotheses.  The  likelihood  ratio  test  for  SIRPs  takes  on  the  form 


w  h2JN(Ql)>  rp 

A(x)  =  w£)  a 


(1.3) 


where  q0  —  xH£  Tx,  qi  =  xH£  *x  -  and 

h2jN(q)=  f  v~2JNexp(—^)fv{v)dv. 
JO  v 


(1.4) 


Using  p2  =  the  magnitude  squared  coherence  between  s  and  x,  q\  is 

expressed  as  qi  =  qo(l  -  (?)■  Since  SIRVs  are  closed  under  linear  transformations,  ro¬ 
tational  and  scaling  operations  on  every  SIRV  result  in  another  SIRV  having  the  same 
characteristic  PDF  [16].  In  all  detection  analyses  involving  colored  noise  spectra,  there 
is  no  performance  penalty  for  whitening  the  noise.  Accordingly,  we  use  an  eigenvalue  de¬ 
composition  of  the  covariance  matrix  of  the  form  £  =  EDEH,  where  E  =  Matrix  whose 
columns  are  orthonormal  eigenvectors  of  £  and  D  is  the  diagonal  matrix  of  eigenvalues 
of  £.  Let  D~V2Effx  =  x'  and  D'^E^s  =  s'.  From  the  closure  property  of  SIRVs  it 
follows  that  x'  is  an  SIRV  with  zero  mean,  identity  (JN  x  JN)  covariance  matrix  and 
characteristic  PDF  /y(u). 

Consequently, 
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where  ||.||2  denotes  the  squared  norm.  , 

We  then  define  an  additional  transformation  A  which  rotates  s  and  places  the  signal 
energy  into  the  first  component  of  the  rotated  vector.  More  precisely  let  s0  =  As  where 
So  =  ||s'||[l,0,...0]r.  The  matrix  A  is  specified  by  the  Householder  transformation  as 
A  =  I  -  2fr£  where  u  =  (s'  -  s0).  Let  x0  =  Ax'.  The  following  properties  of  A  can  be 
readily  established.  For  the  sake  of  completeness,  the  proofs  are  provided  here  although 
they  may  be  found  elsewhere  [17]. 

1.  A  is  a  Hermitian  symmetric  matrix,  i.e.,  A  =  AH. 

Proof: 

A«  =  [I-2{^]»=I-2|f|I=A 

2.  A  is  a  Unitary  matrix.  Thus,  AAH  =  AHA  —  I. 

Proof'  uu^ITt  rjuu^l 

From  property  1,  we  have  A  -  AH.  Hence,  AA  =  A  A  =  [I  — 2pipJ[I  — 2p|pJ  — 


T  / 1  UUg  I  A  UUH  —  T 
1  "  4llulP  +4||uip  ~  1 


3.  The  Householder  transformation  is  norm  preserving. 

We  have  s„  =  As'.  Hence,  ||s0||2  =  s'»A»As'  =  ||s'||2  (since  A»A  =  I  from 
property  2). 

4.  The  inner  product  between  two  vectors  undergoing  the  same  Householder  transfor¬ 
mation  remains  unchanged. 

Proof:  TT 

shXq  =  [As']hAx'  =  s'HAHAx  =  sHx'  (since  AHA  =  I  from  property  2). 

We  use  properties  3  and  4  of  the  Householder  transformation  to  rewrite  q0  and  qx  as 

go  =  ||x'||2  =  ||x0||2  =  E^|Aoi|2  (16) 

9,  =  ||x’||2(i  -  p'2)  =  IM|2(1  -  H.)  =  ES  i*«i2- 

It  is  also  worthwhile  noting  that  x0  is  also  an  SIRV  with  zero  mean,  identity  covariance 
matrix  and  characteristic  PDF  fv{v)-  Let  W  =  |Xoi|2  and  A  —  2  \Xoi\  •  As  a  resu  t, 

the  test  of  eq  (1.3)  reduces  to 


H 1 
> 


A(x)  = 


^2  Jv(A) 
h2jN(W  +  A)^ 


(1.7) 


which  can  equivalently  be  expressed  as 

W>  h^N[^h2JN(A)}  -  A  (1-8) 

Let  T*  =  ho") ar  ( A)]  —  A.  This  reveals  that  the  optimal  strategy  for  detecting  known 

signals  (with  unknown  complex  amplitude)  in  SIRPs  (with  known  covariance  matrix)  is 
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equivalent  to  a  matched  filter  compared  to  a  data  dependent  threshold.  The  test  of  the 
form  of  eq  (1.8)  is  convenient  for  carrying  out  performance  analyses  since  it  enables  us  to 
exploit  the  statistical  independence  between  W  and  A  conditioned  on  V.  This  becomes 
clear  in  the  next  section. 

The  advantage  of  the  Householder  transformation  is  that  it  affords  a  dimensionality  re¬ 
duction  by  using  a  signal  vector  representation,  in  which  all  the  signal  energy  is  contained 
in  the  first  component  of  the  transformed  signal  vector.  The  remaining  components  of 
the  transformed  signal  vector  are  zero.  The  Householder  transformation  also  enables  us 
to  demonstrate  the  maximal  invariance  of  the  test  statistic  of  eq  (1.3)  with  respect  to 
rotation  and  scaling  transformations  and  hence,  establish  the  constant  false  alarm  rate 
(CFAR)  feature  of  the  test. 


1.3  Performance  Evaluation 


In  this  section,  we  derive  integral  expressions  for  the  probability  of  detection  and  false 
alarm  for  the  test  of  eq  (1.8).  Since  x0  is  an  SIRV  with  zero  mean  and  identity  covariance 
matrix,  it  can  be  represented  as  the  product  of  a  complex-Gaussian  random  vector  z  hav¬ 
ing  zero  mean  and  identity  covariance  matrix  and  a  statistically  independent  nonnegative 
random  variable  V  with  PDF  fv(v)  [18].  Thus,  x0  =  zV.  Consequently,  under  the  H0 
hypothesis,  conditioned  on  V,  W  and  A  are  statistically  independent  Chi-Squared  dis¬ 
tributed  random  variables  with  1  and  JN  —  1  complex  degrees  of  freedom,  respectively. 
Furthermore,  the  PDF  of  A  remains  unchanged  under  the  Hi  hypothesis.  Consequently, 

fw\v(w\v)  =  ^exp(-^)  w>0  (1.9) 


A  JN-2  A 

/a|v(A|u)  =  r(JiV  -  1  )u2JV-2eX^('"^2)  A  >  0.  (1.10) 

The  false  alarm  probability,  P/0,  of  the  test  of  eq  (1.8)  conditioned  on  A  and  V  is  given 

Pfa\A,v  =  J™  ^ex P(~)dw  =  exp(~).  (1-H) 

The  unconditional  false  alarm  probability  is  obtained  as 

Pfa  —  Jq  exp(-^)fA\v(&\v)fv(v)dAdv.  (1-12) 

Under  the  Hi  hypothesis,  the  PDF  of  A  conditioned  on  V  remains  unchanged.  However, 
under  Hi,  W  is  recognized  to  be  the  squared  magnitude  of  a  complex  Gaussian  random 
variable  with  mean  rj  =  |a|v'sffE-1s.  Hence,  W  conditioned  on  V  is  simply  the  square  of 
a  Rician  distributed  random  variable  R ,  whose  PDF  is  given  by 


2r 

/«(r)  =  Z2eXPl- 


•£A,(2?)  r  >  0 


(1.13) 
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where  /„(.)  is  the  modified  Bessel  function  of  the  first  kind  of  order  zero.  Hence,  the 
conditional  probability  of  detection  is  given  by 

Pd\a,v  =  J^-fR(r)dr.  (I-14) 

The  unconditional  probability  of  detection  is  obtained  as 

pD=  f°°  I™  PD\A,vfA\v{&\v)fv(v)dAdv.  (1-15) 

Unfortunately,  it  is  not  possible  to  further  simplify  the  integrals  in  eqs  (1.12)  and  (1.15). 
Consequently,  the  integrals  must  be  evaluated  using  numerical  techniques.  We  consider 
the  special  case  of  Gaussian  clutter.  Then,  h2jN{q)  =  exp(—q),  and  fv(v)  =  8(v  -  1). 
Consequently,  T*  =  !£?„[$  WA)1  -  A  =  ln(T)  =  Ta  (say).  For  this  special  case,  the 
data  dependence  of  T*  disappears.  Using  these  observations  in  eq  (1.12)  and  simplifying, 
it  follows  that  Pfa  =  exp(-Ta).  Using  h2jN{-)  =  exp(-.),  fv(v)  =  8(v  -  1)  and  T  =Tq 
in  eq  (1.15),  we  see  that  the  resulting  probability  of  detection  is  given  by 


Pd=  f°°  2rexp[-(r2  +  n2)]I0(2r]r)dr. 

Jy/T * 


(1.16) 


This  is  the  well  known  expression  for  the  matched  filter  detection  probability. 

We  now  generalize  the  analysis  to  allow  for  the  complex  signal  amplitude  to  be  a  random 
variable.  More  precisely,  we  assume  that  a  is  a  complex-Gaussian  random  variable  with 
zero  mean  and  variance  a2.  We  then  derive  the  performance  of  the  test  of  eq  (1.8)  for 
this  case. 

Proceeding  as  before,  it  follows  that  the  PDFs  of  W  and  A  under  the  H0  condition 
remain  unchanged.  Furthermore,  the  invariance  of  the  PDF  of  A  under  Hi  is  also  pre¬ 
served.  It  is  worthwhile  noting  that  the  complex  Gaussian  distribution  on  the  signal 
amplitude  is  also  preserved  after  the  Householder  transformation.  The  only  difference 
under  the  Hi  hypothesis  is  the  conditional  PDF  of  W .  More  precisely,  conditioned  on  V, 
W  is  the  squared  magnitude  of  a  complex-Gaussian  random  variable  with  zero  mean  and 
variance  (5  =  a2sH E^s  +  v2.  Hence,  the  conditional  probability  of  detection  is  given  by 
pD\A  v  =  jexp(-%)dr  =  exp(-y).  The  unconditional  probability  of  detection  is  then 

obtained  by  using  this  expression  in  eq  (1.15).  Again,  as  in  the  previous  case,  closed-form 
evaluation  of  performance  is  not  possible.  The  resulting  integral  must  be  calculated  by 
using  numerical  methods.  The  Gaussian  clutter  case  arises  when  h2JN(x )  =  exp(-x).JLn 
this  case,  fv(v )  =  6{v  -  1).  Using  these  facts  in  eq  (1.15),  it  follows  that  PD  =  exp{—j) 
where  f3  =  a2sH^-ls+l  and  T*  =  ln{T).  As  before,  the  data  dependence  of  the  threshold 

disappears.  .  T  _  , .  ,  , ,  ,  , 

Performance  analysis  is  presented  for  the  case  of  the  K-distribution.  The  K-distnbuted 
envelope  PDF,  which  is  commonly  used  for  modelling  the  amplitude  statistics  of  land  and 
sea  clutter  [19,20],  is  a  specific  case  of  SIRPs.  The  in-phase  and  quadrature  components 
of  a  K-distributed  SIRV  follow  the  generalized  Laplace  distribution.  For  this  case,  the 
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Figure  1.1:  Probability  of  detection  vs  Signal-to-noise-ratio  P/0  =  10  4 

characteristic  PDF  and  h2JN(.)  are  given  by 

fv(v)  =  (6u)2"-1ea;p(-62u2)u(u) 

h2JN(q)  =  2f&-(bVq)a-JNKa-jN(2b^q) 

where  b  is  the  scale  parameter,  a  is  the  shape  parameter,  and  u(.)  is  the  unit  step  function. 
We  use  fv(v)  and  h2JN{.)  for  the  K-distribution  in  eq  (1.12)  to  determine  the  threshold 
for  a  false  alarm  probability  of  10-4.  Detection  probability  is  calculated  using  eq  (1.15) 
by  numerical  integration.  Probability  of  detection  is  plotted  as  a  function  of  the  signal- 
to-noise-ratio,  |d|2sHS_1s,  for  the  case  of  a:  =  0.5  and  a  =  0.1  in  Figures  1.1  and  1.2 
respectively.  Relevant  test  parameters  are  reported  in  the  figures.  In  [21],  the  performance 
evaluation  of  the  optimal  receiver  was  carried  out  via  Monte-Carlo  simulations.  The 
results  of  our  analysis  for  the  two  test  cases  considered  are  in  agreement  with  the  results 
of  [21]. 


1.4  Conclusions 

We  considered  the  performance  of  the  optimal  processor  for  target  detection  in  additive 
non-Gaussian  clutter  which  can  be  modeled  as  an  SIRP.  We  derived  integral  expressions 
for  the  probability  of  false  alarm  and  probability  of  detection  for  the  case  of  a  known 
signal  with  unknown  complex  amplitude.  The  results  were  then  extended  to  the  case  of  a 
signal  with  unknown  complex  amplitude  following  a  complex-Gaussian  distribution.  The 
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Figure  1.2:  Probability  of  detection  vs  Signal-to-noise-ratio  Pfa  -  10  4 

Householder  transformation  provided  a  useful  technique  for  simplifying  the  analysis  As 
a  result,  several  important  properties  of  the  Householder  transformation  have  been  doc- 
umented  in  this  work.  Our  results  reduce  to  the  well  known  expressions  for  the  Gaussian 
matched  filter  as  a  special  case.  Except  for  the  special  case  of  Gaussian  clutter,  closed 
form  or  near  closed  form  evaluation  of  the  integrals  is  not  possible.  Further  work  needs 
to  be  undertaken  to  improve  the  numerical  integration  method  used  in  this  work. 
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Chapter  2 

Analytical  Expressions  for  the  NMF 
and  NAMF  Tests 

2.1  Test  Statistic  Descriptions 

We  now  consider  several  non-adaptive  and  adaptive  detection  test  statistics  in  this  section. 


2.1.1  Non- Adaptive  Test  Statistics 

The  optimal  receiver  for  detecting  a  rank  one  signal  in  Gaussian  interference  with  known 
covariance  matrix,  R is  given  by 


A  mf  — 


le^R^xl2 

e^R^e 


Hi 

<  A  mf- 

Ho 


(2.1) 


In  some  instances,  the  test  data  vector  is  subject  to  an  unknown  scaling,  if.  Consequently, 
the  test  data  covariance  matrix  takes  the  form  fUd-  This  corresponds  to  the  scenario, 
where  the  covariance  matrix  has  known  structure  but  unknown  level.  The  phase  invariant 
matched  filter  (PI-MF)  test  for  these  problems  is  expressed  as  [22] 


A  PIMF  — 


|e*R^x|2  g 
rj2eH Rj*e  h0 


A PIMF 


(2.2) 


where  e  and  x  are  the  concatenated  JN  x  1  signal  ‘search’  steering  and  data^  vectors, 

respectively.  The  inner  product  of  whitened  vectors  b  =  R^^x  and  f  =  R<;  is  the 
matched  filtering  operation.  Although  (2.2)  does  not  require  knowledge  of  signal  phase,  it 
does  require  knowledge  of  the  scale  77  to  be  CFAR.  The  normalized  matched  filter  (NMF) 
test  [14, 22-28]  is  given  by 


AatmF  — 


|e*Rd-] 


[etfR/ejfxtfR/x] 


Hi 

<  Aatmf- 

Ho 


(2.3) 
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Th-  whitening  operation  R?  transforms  the  signal  and  data  vectors  to  obtain  f and lb, 

Schwartz  inequality,  0  <  A mf  < 

2.1.2  Adaptive  Test  Statistics 

For  the  adaptive  problem,  Rs  replaces  R„  and  the  tests  of  (2.1)  and  (2.3)  are  referred  to 
here  as  the  AMF  and  NAMF,  respectively.  These  tests  are  given  y 


A  amf  — 


|C  >  XNAMF. 

Wr  -  letfRj'ellx^'x]  i. 

Observe  that  A  amf  is  simply  the  adaptive  version  of  A„f  for ;  ^ 

«  hnown  as 

the  Kelly  GLRT  [6]  is  expressed  as 

le^xl2  H>Vi_  12.6) 


i^‘r2>w 

(2.4) 

eHRJ1e  »o 

>' WF. 

(2.5) 

*  „  _  >K\GLkt  W 

KaLRT  ~  .  ,  „  x"R,71x,  «. 

[e«R;'e][l  +  —^—1 

where  0  <  AGiKr  <  1.  It  la  instructive  to  note  that  in  the  limit  of  large  K,  the  tests  of 
!£££ ZC^tZ  in  ^Jimit  XXhe  tests  Converges  to 

noted  to  be  CFAR  tor  the  prooiem  oi  «  2  .  while  the  training  data  z 

No  optimality  claims  of  the  NAMF  test  can  be  made  for  the  case  of  SIRP  disturbance. 

2.2  Analytic  Results 

-sr,r^ 

considered  for  tractability  of  analysis. 
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The  test  statistic  of  (2.3)  can  be  expressed  in  terms  of  the  whitened  steering  vector,  f, 
and  the  whitened  test  data  vector,  b  as 


A2 


|fKb|2  "1 
“  [ffff][bHb] 


A2 


(2.7) 


Noting  that  a  unit  vector  in  the  direction  of  f  is  given  by  fi  =  (f4)o.s>  bKb  can  be 
expressed  as  the  sum  of  the  squared  magnitudes  of  projections  along  the  subspace  of  f 
and  the  orthogonal  complement  space  of  f  denoted  by  P_l-  Let  w*,  i  =  1,2, ...  JN  —  1 
denote  an  orthonormal  basis  set,  for  Pj_  and  Xq  =  f^b,  Xi  —  w^b,  i  =  1,2, ...  JN  —  1. 
Then,  Xu  i  =  0, 1, . . .  JN  -  1  are  mutually  uncorrelated  univariate  complex-SIRVs.  Let 

JN- 1 

&  =  |x0|2,  6  =  S  l^l2’  and  $  =  f;-  The  test  statistic  of  (2.7)  admits  a  representation 

Z=1 


of  the  form 


2  (1  +  $) 


(2.8) 


Under  H0,  Xi,  i  =  0, 1, . . . ,  JN  -  1  are  zero  mean  unit  variance  complex-SIRVs  with 
characteristic  PDF  fv(v )•  Hence,  Xi  =  ZJJ ,  where  Zi,  i  —  0, 1, . . . ,  JN  —  1  ,  are  iid 
CN(0,1)  random  variables.  Let  Xi  =  \%o\2 ■>  Xjn-  1  =  J2i=i~l  \Zi\2  >  and  A  =  —  . 

XJN—1 

Thus,  under  H0,  $  =  A.  The  test  statistic  of  (2.7)  reduces  to 


A2  — 


A 

(1  +  A) 


(2.9) 


A  is  the  ratio  of  two  statistically  independent  Chi-Squared  distributed  random  variables. 
Hence  it  follows  the  central-F  distribution  [30]  given  by 


•k®  0(JN  -  1, 1)(1  +  5)JN 


(2.10) 


where  x 

/ 3(m,n)=  [  (2.11) 

JO 

From  (2.9)  using  a  straightforward  transformation  of  random  variables,  it  follows  that  the 
PDF  of  A2(x)  under  Hq  is  given  by 

hJr)  =  (JN-l)(l-r)JK-2.  (2.12) 


Consequently,  the  probability  of  false  alarm  is  given  by 

Pfa-NMF  =  R(A2  >  A2|i7o)  =  (1  —  A 2)JN  X-  (2.13) 

Under  Hu  conditioned  on  V,  &  is  a  non-central  Chi-Square  distributed  random  variable 
with  noncentrality  parameter  A  =  |a|(e^R^’1e)  2  while  £2  conditioned  on  V  is  a  central  chi- 
square  distributed  random  variable  with  JN  —  1  degrees  of  freedom.  Hence,  conditioned 
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on  V,  $  has  a  non-central  F  distribution.  After  some  algebra,  the  PDF  of  A2  given  V  is 
expressed  as 

/  ^  A2krk{l-r)JN-2  (0-iA\ 

f\2\v(r\v)  =  T,exp(-^  )~2i^(jN  -  1,  kTT)  °-r-  (‘  } 

The  conditional  probability  of  detection  is  given  by 

00  ^^2  ^2/c  \*l 

Pi-HMFV  =  F[( A2|V)  >  A2|ff,]  =  -  lxtainc(\2,k  +  1,JN-1)} 

k~°  (2.15) 

where  -1  rX  , 

>>etainc(x,  m,  n )  =  ^^y/0  C”'1*1  -  <)”'«•  <2'16) 

The  unconditional  probability  of  detection  obtained  by  taking  the  expectation  of  (2.15) 

over  V  is  given  by 


oo  A2k 

Pd-NMF  =  -  tetainc(\2,  k  +  1,  JN  -  1)] 


(2.17) 


where 


h2k(<P)  =  f  «  2texp(-^)fv(v)dv. 


(2.18) 


Finite  sum  expressions  for  detection  probabilities  are  reported  in  [31].  Using  the  results 
therein,  the  expression  of  (2.17)  simplifies  to 


r  {JN) 


Pd-NMF  - 1  -  (i  -  a2 )jn  1  E  +  i)r (jn  -k){ i  -  a2 


where 


Wl 


gk[A\\  -  A2)]  (2.19) 


9‘W  =  §r^T)M'")-  <2'20> 

For  the  case  of  the  K-distributed  SIRP, 

9  fc-1  /3/+a?/.0.5a+0.5i  _ 

=  rRS^WTTT^-''2^  (2'21) 

For  convenience,  we  set  0  =  y/a.  in  the  examples  presented.  This  normalization  ensures 
unit  mean  square  value  for  the  modulating  random  variable  V.  The  result  of  (2.19)  is 
used  in  the  examples  presented  in  the  next  section.  For  the  special  case  where  h2l{w)  - 
exp(—w) , (corresponding  to  Gaussian  clutter),  the  expression  of  (2.19)  reduces  to 


r  {JN)  f  A2_V 

Pd-NMF  =  l— (i— a2)  E  +  i)r(jAr  -k){ i  -  a2 J 


[1— gammainc(A2  (1— A2),  k )] 
(2.22) 


11 


where 


1  cp 

gammainc(p,m )  =  p^yyo  0)d6 


(2.23) 


is  the  incomplete  Gamma  function. 

We  now  consider  the  problem  of  performance  evaluation  of  the  NAMF  in  Gaussian 
clutter.  The  NAMF  test  statistic,  A NAmf  is  related  to  the  Kelly  GLRT  test  statistic, 

A rr  jrr  of  (2  6)  where  T  = - l— - is  the  well  known  loss-factor  with  PDF  [5] 

v  '  [l  +  x#R J  x/K] 


/r(7>  =  /?(L  + 1,  JJV -T)7^1  ~  'r)JN~2 


(2.24) 


where  L  =  K  -  JN  +  1.  Noting  that  Anamf  =  f°U°WS  fr°m  ^  that  the 

probability  of  false  alarm  conditioned  on  F  is  given  by 


P fa— NAM  F I F 


[l  +  (l  -  'y)v]i 


(2.25) 


where  rj  is  the  threshold.  The  unconditional  probability  of  false  alarm  is  obtained  by 
taking  the  expectation  of  (2.25)  over  T  and  is  given  by 


p'-“"  =  L  [i  +  (i-7 wdr  1 

Similarly  from  [6]  the  NAMF  probability  of  detection  conditioned  on  T  is  given  by 
-  l-nTTI^grt  (  i  )  m-lTHsammainc 


Mi) 


(2.26) 


[rj  ( 1 7)  ] m  [1  —  gammainc 


=  1-[1  +  (1_T)#Z  (  m  )  W-lT^Samma inc  (j -  ^™)l- 

l  m-1  (2.27) 

Taking  the  expectation  of  (2.27)  over  F  yields  the  unconditional  probability  of  detection, 
which  is  expressed  as 

=  1-/„1[1  +  (11_t)#£  (  m  )  [»(l-T)ni-9— no 

(2.28) 

For  convenience,  Pd-NAMF  is  expressed  in  terms  of  P<icfar  of  a  scalar  CFAR  processor 
given  by 


-NAMF 


=  i-/1 _ I - V) 

io  [l  +  (l-7)#ii 


Pd  cfar(v,  A,L)  =  1-  ^2--  E  ( t )  ^  -  »ammainc  (tT7  m)  1-  (2  29) 


Numerical  calculation  of  Pd-NAMF  is  obtained  as 


M-1 

Pd-NAMF  =  A7  PdCFax[(l  —  ^7)77,  (1  —  ^ A7)A,  T]/r(£ A7) 


(2.30) 
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where  A7  =  4j  .  Analytic  expressions  for  the  Pfa  and  Pd  of  the  NAMF  operating  in 
Gaussian  clutter  were  independently  derived  in  [2,32].  The  corresponding  expressions  for 
the  NMF  in  compound-Gaussian  clutter  appear  here  for  the  first  time. 

In  general,  it  is  extremely  difficult  to  derive  the  corresponding  Pfa  and  Pd  for  the 
NAMF  operating  in  SIRP  clutter.  This  is  due  to  the  fact  that  the  ML  estimate  of  the 
SIRP  covariance  matrix  is  not  available  in  closed  form  [33].  Hence,  its  PDF  cannot 
be  determined  analytically.  For  the  special  case  where  the  ML  estimate  of  the  SIRP 
clutter  covariance  matrix  is  to  within  a  multiplicative  constant  of  the  sample  covariance 
matrix  [34],  it  is  possible  to  derive  the  P/a  and  Pd  expressions  for  the  NAMF  in  SIRP 
clutter.  Specifically,  for  this  scenario,  the  expression  for  Pfa  remains  unchanged  from 
(2.26).  However,  the  result  for  Pd  is  quite  different  from  that  of  (2.28).  The  scenario 
where  the  ML  estimate  is  to  within  a  multiplicative  constant  of  the  sample  covariance 
matrix  corresponds  to  the  case  where  the  texture  component  of  the  SIRP  is  perfectly 
correlated.  More  precisely,  for  this  scenario  the  texture  component  remains  constant  over 
all  the  training  data  realizations.  However,  real  data  seldom  exhibits  such  a  behavior  [35]. 
Consequently,  Pfa  and  Pd  results  for  this  scenario  are  merely  of  academic  interest  having 
little  or  no  practical  importance. 

The  Pfa  and  Pd  expressions  for  the  N-PAMF  and  PAMF  operating  in  both  Gaussian 
and  non-Gaussian  clutter  scenarios  are  lacking  since  the  analysis  becomes  mathematically 
intractable.  Additionally,  it  is  difficult  to  guarantee  the  CFAR  behavior  of  the  N-PAMF 
and  PAMF  with  respect  to  the  unknown  covariance  matrix.  Consequently,  performance 
evaluation  of  these  methods  is  carried  out  by  Monte-Carlo  techniques.  Peformance  results 
are  presented  in  the  next  section. 


2.3  Performance  Results 

Performance  is  now  presented  for  the  detectors  described  above.  Figure  2.1  depicts  the 
performance  of  the  NMF  in  terms  of  Pd  versus  SINR  in  K-distributed  SIRP  for  several 
shape  parameter  values.  Relevant  test  parameters  are  provided  in  the  plot.  The  Pd  curves 
shown  in  this  figure  are  obtained  by  using  (2.19)  and  (2.22).  It  is  instructive  to  observe  the 
potential  for  improved  NMF  performance  in  impulsive  clutter  (a  =  0.1).  Furthermore, 
the  NMF  performance  represents  the  upper  bound  on  the  performance  of  the  N-PAMF 
and  NAMF. 

Figure  2.2  plots  the  NMF  threshold  as  a  function  of  the  clutter  shape  parameter. 
Observe  that  the  threshold  is  independent  of  the  clutter  shape  parameter.  Hence,  the 
NMF  is  texture  CFAR.  This  property  is  of  considerable  importance  in  practice,  where 
estimation  of  the  clutter  shape  parameter  imposes  onerous  requirements  of  training  data 
support.  However,  this  property  of  the  NMF  is  lost  when  the  covariance  matrix  underlying 
the  clutter  is  unknown.  Specifically,  it  is  shown  in  [36]  that  while  operating  SIRP  clutter 
the  NAMF  employing  a  sample  covariance  matrix  estimate  incurs  a  considerable  increase 
in  the  false  alarm  probability  for  small  values  of  the  shape  parameter. 

Figure  2.3  shows  the  performance  of  the  NAMF  operating  in  Gaussian  clutter.  Also 
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Figure  2.3:  NAMF  Performance  in  Gaussian  Clutter 


shown  is  the  performance  of  the  NMF  in  Gaussian  clutter.  We  observe  that  with  K— 2JN, 
for  a  given  probability  of  detection,  the  SNR  required  by  the  NAMF  is  approximately  3 
dB  larger  than  the  corresponding  SNR  for  the  NMF . 

Figure  2.4  presents  a  performance  comparison  of  the  NAMF,  AMF  and  Kelly  GLRT 
in  Gaussian  clutter.  For  the  example  presented,  the  Kelly  GLRT  performs  slightly  better 
than  the  NAMF  and  CFAR-AMF.  The  latter  two  methods  have  comparable  performance. 


2.4  Summary 

In  this  chapter,  we  presented  analytic  expressions  for  the  probability  of  detection  and 
probability  of  false  alarm  for  the  NMF  operating  in  compound-Gaussian  clutter.  The 
texture  CFAR  feature  of  the  NMF  was  noted.  Probability  of  detection  and  false  alarm 
for  the  NAMF  operating  in  Gaussian  clutter  is  also  presented.  Performance  of  the  NAMF 
is  compared  with  that  of  the  NMF  operating  in  Gaussian  clutter.  Approximately  3dB 
difference  in  performance  was  noted  for  K=2JN. 
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Probability  of  detection  (Pd) 


Chapter  3 


Statistical  Analysis  of  the 
Nonhomogeneity  Detector 

3.1  Introduction 

An  important  issue  in  space-time  adaptive  processing  (STAP)  for  radar  target  detection  is 
the  formation  and  inversion  of  the  covariance  matrix  underlying  the  clutter/interference. 
In  practice,  the  unknown  interference  covariance  matrix  is  estimated  from  a  set  of  in¬ 
dependent  identically  distributed  (iid)  target-free  training  data  which  is  assumed  to  be 
representative  of  the  interference  statistics  in  a  cell  under  test.  Frequently,  the  training 
data  is  subject  to  contamination  by  discrete  scatterers  or  interfering  targets.  In  either 
event,  the  training  data  becomes  nonhomogeneous.  As  a  result,  it  is  non  representative  of 
the  interference  in  the  test  cell.  Estimates  of  the  covariance  matrix  from  nonhomogeneous 
training  data  result  in  severely  undernulled  clutter.  Consequently,  CFAR  and  detection 
performance  suffer.  Significant  performance  improvement  can  be  achieved  by  employing 
pre-processing  to  select  representative  training  data. 

The  problem  of  target  detection  using  improved  training  strategies  has  been  considered 
in  [37-44].  The  impact  of  nonhomogeneity  on  STAP  performance  is  considered  in  [44-47]. 
The  works  of  [37-40,44,48]  have  addressed  the  use  of  the  non-homogeneity  detector  (NHD) 
based  on  the  generalized  inner  product  (GIP)  measure  for  STAP  problems.  The  GIP  is 
used  as  a  measure  to  select  representative  training  data  by  the  NHD.  In  this  chapter,  we 
provide  a  statistical  analysis  of  the  GIP  based  NHD  for  Gaussian  interference  statistics. 


3.2  GIP  Statistics:  Known  Covariance 

Let  x  =  [xi  X2  ...  xm]t  denote  a  complex  random  vector  with  zero  mean  and  known 
positive  definite  Hermitian  covariance  matrix  R.  The  quadratic  form  given  by  Q  — 
xHR_1x  has  the  important  property  that  E(Q)  =  M  [30].  This  result  is  readily  proven 
below. 

Q  =  xhR_1x  =  fr[xffR-1x]  =  tr[R_1xxH]  (3.1) 
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where  tr(.)  denotes 
tr(BA).  Hence, 


the  trace  of  a  matrix.  We  have  made  use  of  the  fact  that  tr( AB)  — 

E(Q)  =  E{tr[ R^xx*]}  =  frfRr^xx*}]  /3  2) 

=  tr(  RXR)  =  tr(IM)  =  M 


where  IM  is  the  M  x  M  identity  matrix.  This  result  is  important  in  that  it  is  independent 
of  the  PDF  underlying  x  and  is  only  a  function  of  the  dimension  of  the  random  vector. 

If  the  PDF  of  x  is  known,  the  corresponding  PDF  of  Q  can  be  readily  derived.  For 
Gaussian  distributed  x,  i.e.,  x  ~  CN( 0,R),  the  PDF  of  Q  is  a  Chi-Squared  distribution 
with  M  complex  degrees  of  freedom.  More  precisely, 


Q  =  xhR-1x  =  ||R  2x||2 


(3.3) 


where  ||.||  denotes  the  Euclidean  vector  norm.  Letting  y  —  R  5x  gives 

M 

<?=l|y||2  =  £WI2 


(3.4) 


where  Yi,  i  =  1, 2, . . .  M  are  iid  CN( 0, 1)  random  variables.  Since  Q  is  the  sum  of  the 
squared  magnitudes  of  M  iid  CN{ 0, 1)  random  variables,  it  follows  that  Q  is  a  Chi 
Squared  distributed  random  variable  with  M  complex  degrees  of  freedom  [30].  The  PDF 
of  Q  is  given  by 

n  M— 1  v 

/q(?)  =  f{M)eXp ^  °-q<0° 

where  T(.)  is  the  Eulero-Gamma  function. 

The  GIP  based  NHD  calculates  the  quadratic  form  Q  using  an  estimated  covariance 
matrix  (formed  from  iid  target  free  training  data)  and  compares  its  mean  with  the  di¬ 
mensionality  of  the  random  vector  x.  Deviations  from  M  have  been  attributed  to  non¬ 
homogeneities  in  the  training  data  [37-40,44],  In  practice,  the  interference  covariance 
matrix  is  formed  from  a  finite  amount  of  training  data.  The  statistical  variability  associ¬ 
ated  with  the  data  could  introduce  additional  errors  and  thus  deviations  of  the  GIP  from 
M  cannot  entirely  be  ascribed  to  the  presence  of  non-homogeneities.  Consequently,  it  is 
useful  to  work  with  the  statistics  of  Q  formed  with  an  estimated  covariance  matrix. 


3.3  GIP  Statistics:  Unknown  Covariance 

The  complex-Gaussian  test  data  vector  is  denoted  by  x  ~  CN(0,Ht),  where  Rt  is 
unknown.  Let  Zi,  i  =  1,2, . .  .K  denote  iid  CN( 0,R*)  target  free  training  data.  For 
representative  (homogeneous)  training  data,  RT  =  R *  =  R-  The  maximum  likelihood 
estimate  of  the  covariance  matrix  is  given  by 


1* 


H 


(3.6) 
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p 

Let  P  =  xFR_1x.  A  stochastic  representation  for  the  normalized  GIP  denoted  by  P  =  — 
is  derived  in  appendix  A.  Consequently,  P  can  be  expressed  as 


P  = 


KRi 

P‘2 


(3.7) 


where  Ri  and  R2  are  statistically  independent  Chi-squared  distributed  random  variables 
with  PDFs  given  by 


rM- 1 

/ft  (rl)  =  jh gyeip(-n)  0  <  r,  <  00 


(3.8) 


K—M 

/ifafa)  =  r  (KlM  +  t)eXp^r^  °-r2<°°  (3-9) 

respectively.  The  PDF  of  P' ,  which  is  simply  a  central-F  distribution  [49],  is  expressed  as 

1  rM-l 

fp'  ^  “  0(L,  M)  (1  +  r)M+L  (3.10) 

0  <  r  <  oo 


where  L  =  K  -  M  +  1  and 

p(m,n )  =  /V-^l  -  e)n~ldO.  (3.11) 

The  statistical  equivalence  of  P  to  within  a  scalar  of  the  ratio  of  two  independent  chi- 
square  distributed  random  variables  is  fascinating  in  that  it  permits  rapid  calculation 
of  the  moments  of  P.  More  importantly,  it  is  extremely  useful  in  Monte-Carlo  studies 
involving  computer  generation  of  P.  For  homogeneous  training  data,  the  use  of  (3.7) 
circumvents  the  need  to  explicitly  generate  the  test  data  vector  x  and  the  training  data 
vectors  used  for  covariance  estimation.  For  large  M  and  perforce  K,  significant  compu¬ 
tational  savings  can  be  realized  from  the  method  of  (3.7). 

It  can  be  readily  shown  that 


E(P)  =  KE(R1)E(R21) 


-i\  - 


M 


4  =  7 


M 


x-E 

K 


(M±l) 

1  K  1 


(3.12) 


where  E(P)  and  aj>  denote  the  mean  and  variance  of  P,  respectively.  We  then  consider 

the  PDF  of  —  in  the  limit  of  large  K.  The  characteristic  function  of  —  is  given  by 
K  ^ 

$n2{juj)  =  E[exp(-ju^)}  = - jJ- - •  (3.13) 

K  (1  +  ^)^"M+1 
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(3.14) 


In  the  limit  of  K  -4  oo,  we  have 

tf(jw)  =  limK^oo$R2(ju)  =  exp(-ju). 


R 

The  PDF  of  -y  in  the  limit  of  K  -»  oo  is  given  by 
K 

fR2(r )  =  y-  f°°  ercM(r  ~  l)M<^  =  “  !)•  (3-15) 

Hence,  for  K  ->  oo,  tf2/A  becomes  unity  with  probability  1.  Consequently,  the  GIP 
is  simply  Rt  and  hence,  follows  a  Chi-Squared  distribution  with  M  complex  degrees  of 
freedom.  Thus  for  K  -4  oo,  E(P)  =  M  and  a%  =  M,  corresponding  to  the  known 
covariance  matrix  results.  Consequently,  the  GIP  statistical  representation  given  by  (3.7) 
provides  additional  insights  into  the  NHD.  The  numerator  random  variable  corresponds 
to  the  GIP  statistics  for  known  covariance  matrix.  The  denominator  random  variable 
succinctly  embeds  the  deleterious  effects  of  estimating  the  covariance  matrix  with  finite 
sample  support  as  well  as  nonhomogeneity  of  the  training  data.  A  manifestation  of 
this  effect  can  be  seen  from  the  deviation  of  the  statistics  of  i?2  from  the  Chi-Squared 
distribution. 


3.4  New  Test  for  Nonhomogeneity 

The  work  of  [37-40,44]  uses  an  NHD  based  on  comparing  the  mean  of  empirically  formed 
GIPs  (from  different  realizations  of  test  data)  to  M.  Large  deviations  from  M  are  as¬ 
cribed  to  nonhomogeneities.  However,  the  effects  of  finite  data  support  and  the  associated 
statistical  variability  can  result  in  large  deviations  of  the  empirical  GIP  mean  from  M. 
Thus,  a  more  stringent  test  for  the  GIP  based  NHD  consists  of  the  following  steps: 

1.  Form  the  GIP  denoted  by  P  for  a  given  K. 

2.  Perform  a  goodness-of-fit  test  of  the  empirically  formed  f  with  the  theoretically 
predicted  PDF  of  (3.10).  Specifically,  we  set  the  type-I  error,  a,  to  be  0.1.  This  is 
simply  the  probability  of  incorrectly  rejecting  the  hypothesis  that  the  data  comes 
from  the  F-distribution  of  (3.10).  More  precisely,  this  corresponds  to  calculating  a 
threshold  A,  such  that  a  =  Pr(P'  >  A)  =  0.1.  From  (3.10),  it  follows  that 

Pr(P'  >  A)  =  betainc(  — ,  M,  L)  (3.16) 

A  +  1 

where  ,  x 

betainc(x,  M,  L)  =  ^ ^  ~^)L  ^ *  (3-17) 

Given  a,  M,  and  L,  A  is  readily  calculated  from  an  inversion  of  (3.16).  The  goodness- 
of-fit  test  consists  of  comparing  the  empirically  formed  P'  from  each  training  data 
realization  with  A  and  rejecting  those  realizations  for  which  P  exceeds  A. 
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3.  Compare  the  empirical  GIP  mean  to  the  theoretically  predicted  mean  value  of  (3.12). 
Large  discrepancies  between  the  theoretical  and  empirical  means  result  from  non¬ 
homogeneities.  This  serves  as  an  independent  validation  of  step  2. 

4.  Histogram  plots  of  the  empirical  values  of  £  overlaid  on  the  theoretical  PDF  of  P' 
serve  as  an  additional  check  on  the  results  of  steps  2  and  3. 

3.5  Performance  Analysis 

We  present  performance  results  of  our  approach  here.  Figure  3.1  shows  the  PDF  of  P 
for  several  values  of  K  with  Ivl — 8  for  Gaussian  interference  statistics.  ^Observe  that  the 
variance  of  P'  decreases  with  increasing  K.  This  is  anticipated  since  R  tends  to  R  with 
probability  1  as  K  ->  oo.  The  results  presented  in  Figures  3.2  and  3.3  correspond  to  the 
case  of  homogeneous  training  data.  Figure  3.2  presents  a  comparison  of  the  cumulative 
distribution  function  (CDF)  of  P  obtained  from  Monte-Carlo  realizations  using  simulated 
data  with  the  theoretically  predicted  CDF  of  P'  obtained  by  numerical  integration  of 
(3.10).  Figure  3.3  plots  the  theoretical  PDF  of  P'  given  by  (3.10)  over  the  histogram 
obtained  from  Monte-Carlo  data  for  M=8  and  K=16.  The  results  show  good  agreement 
between  the  theoretical  prediction  and  the  empirically  generated  values.  The  mean  value 
of  P,  15.957,  obtained  via  Monte-Carlo  compares  well  with  the  theoretically  predicted 
value  of  16.  Figure  3.4  plots  the  type-I  error  versus  threshold  for  M=64.  Here  different 
values  of  K  are  chosen  to  illustrate  the  threshold  behavior.  For  each  value  of  a ,  A  is 
determined  from  a  numerical  inversion  of  (3.16).  For  a  given  a  we  observe  an  increase 
in  A  with  increasing  K.  Figure  3.5  shows  the  results  of  the  goodness  of  fit  test  for  the 
MCARM  data  [50]  using  acquisition  ‘220’  on  Flight  5,  cycle  ‘e’  for  8  channels  and  16 
pulses.  The  normalized  GIP  and  the  threshold  are  plotted  as  a  function  of  range.  Non¬ 
homogeneity  of  the  training  data  is  evident  in  those  bins  for  which  the  normalized  GIP 
exceeds  the  threshold. 


3.6  Conclusion 

This  work  extends  the  results  of  previous  work  on  non-homogeneity  detection  by  providing 
a  rigorous  statistical  characterization  of  the  NHD  for  Gaussian  interference.  It  is  shown 
that  the  NHD  statistic  admits  a  simple  representation  as  a  ratio  of  two  statistically  in¬ 
dependent  Chi-squared  distributed  random  variables.  A  formal  goodness-of-fit  test  based 
on  this  representation,  which  follows  an  F-distribution,  is  derived.  Performance  analysis 
of  the  method  is  considered  in  some  detail.  The  illustrative  examples  validate  the  ap¬ 
proach  taken  and  confirm  the  results.  Future  work  would  include  extensive  performance 
analysis  using  simulated  and  measured  data  showing  the  resulting  impact  on  STAP  per¬ 
formance.  The  performance  of  several  STAP  algorithms  in  Gaussian  and  non-Gaussian 
interference  scenarios  has  been  considered  in  [36].  Future  work  will  address  performance 
of  the  methods  treated  in  [36]  with  suitable  NHD  pre-processing. 
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1.4 


PDF  of  P  for  Several  Values  of  K 


Figure  3.1:  PDF  of  the  normalized  GIP 


Theoretical  and  Empirical  CDF  of  P 


Number  of  Occurences 
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Appendix  A: Stochastic 
Representation  for  the  Normalized 


GIP 


Let  Z  denote  a  data  matrix  whose  columns  are  the  previously  defined  z;,  i  =  1,2,. 
The  maximum  likelihood  estimate  of  the  covariance  matrix  is  then  expressed  as  R  = 
where  S2  =  ZZH.  Consequently,  the  normalized  GIP  is  expressed  as 


..,K. 

Is 


P 


=  xHs;xx. 


(3.18) 


The  data  matrix  Z  and  the  vector  x  admit  a  statistical  representation  of  the  form 


Z  =  R*Y 

_  1 

x  =  R2y 


(3.19) 


where  Y  is  a  data  matrix  whose  columns,  y,,  i  =  1,2,..., A  are  iid  CN(0,I)  random 
vectors  and  y  is  a  CN(0,I)  random  vector,  which  is  statistically  independent  of  Y.  Hence, 
the  normalized  GIP  is  expressed  as 

P  =  y^y  (3.20) 

uuff 

where  Su  =  YYH .  Next,  we  use  a  Householder  transformation  defined  by  A  =  1-2— 

y  u  u 

where  u  =  y  -  ||y||e  and  e  =  [100 . . .  Of,  so  that  y  =  Ay  =  ||y||e.  Also,  let  Y  =  AY. 
Since  A  =  AH  and  AAH  =  AHA  =  I,  it  follows  that  the  statistics  of  Y  are  identical  to 
that  of  Y.  Conseqently,  the  normalized  GIP  is  expressed  as 


p'  =  yHS-ly 

(3.21) 

where  =  YYH.  Furthermore,  we  partition  Y  as 

i - 1 

1  - 1 

II 

«!* 

(3.22) 

29 


where  is  the  the  first  row  of  Y  and  Y^  denotes  the  ( M  —  1)  x  K  matrix  formed  from 
the  remaining  rows  of  Y.  Consequently,  Sy  is  expressed  as 


c  _  f  yfyi  yfYn 
5  [  Ygyi  YgYu 


(3.23) 


The  inverse  of  Sy  admits  a  representation  of  the  form 


(3.24) 


Finally,  the  normalized  GIP  is  expressed  as 


P  = 


(3.25) 


However,  from  the  matrix  inversion  Lemma  it  follows  that  S11  =  [yfPj.yi]-1,  where 
P±  =  [I  -  Yu(YgYii)_1Yg].  Since  Yi^YgYn)"1^  is  a  projection  matrix  of  rank 
M-l,  it  follows  that  P_l  is  a  projection  matrix  of  rank  K-M+l.  Consequently,  yfP±yi  = 

K-M+l 

Y  |yx(i)|2.  Therefore,  S11  is  simply  the  reciprocal  of  a  chi-squared  distributed  random 

i=l 

variable  with  K-M  +  l  complex  degrees  of  freedom.  Also,  ||y||2  is  simply  the  sum  of 
the  squared  magnitudes  of  M  iid  CN(0,1)  random  variables  and  hence  follows  a  chi- 
squared  distribution  with  M  complex  degrees  of  freedom.  Consequently,  the  GIP  admits 
a  representation  of  the  form  of  (3.7). 


30 


